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We analyze critically the renormalization of quantum fields in cosmological spacetimes, 
using non covariant ultraviolet cutoffs. We compute explicitly the counterterms necessary 
to renormalize the semiclassical Einstein equations, using comoving and physical ultraviolet 
cutoffs. In the first case, the divergences renormalize bare conserved fluids, while in the 
second case it is necessary to break the covariance of the bare theory. We point out that, 
in general, the renormalized equations differ from those obtained with covariant methods, 
even after absorbing the infinities and choosing the renormalized parameters to force the 
consistency of the renormalized theory. We repeat the analysis for the evolution equation 
for the mean value of an interacting scalar field. 
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I. INTRODUCTION 


There is a well defined and rigorous approach for the renormalization of quantum fields in 
curved spacetimes [l|. A covariant regularization of the theory (for example point splitting or 
dimensional regularization) is combined with the Hadamard structure of the two point function 
in order to perform adequate subtractions and obtain renormalized expressions for the effective 
action and the mean value of the Energy Momentum Tensor (EMT), the source in the so called 
Semiclassical Einstein Equations (SEE). The infinities of the theory are absorbed into the bare 
constants of the classical gravitational action (Newton constant G, cosmological constant A, and 
dimensionless coefficients of terms quadratic in the curvature). The whole procedure preserves the 
covariance of the theory. 

For different reasons, there have been attempts to regularize the theory using ultraviolet (UV) 
cutoffs. For example, in basic discussions of the cosmological constant problem, a three dimensional 
cutoff is considered in order to obtain a naive estimation of the zero point energy of quantum 
fields, which turns out to be an enormous contribution, and would require a fine tuning of the bare 
cosmological constant of more than 120 orders of magnitude Q]. 

The use of UV cutoffs is more intuitive than dimensional regularization, but has been disre¬ 
garded in semiclassical gravity because of its clash with the covariance of the theory Q]. How¬ 
ever, more recently it has been argued that it is possible, in principle, to renormalize the theory 
in Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetimes using a three dimensional physical 
cutoff 1 ^. Being proportional to the scale factor, the time dependence of the physical cutoff spoils 
the conservation of the EMT and, in order to restore it, it would be necessary to introduce non 
covariant counterterms, whose finite parts should be carefully chosen to maintain the consistency 
of the renormalized theory. Some relevant aspects of the proposal have not been worked out, like 
the calculation of the counterterms for general metrics, and a comparison with the usual approach. 
Regarding the quartic divergences that dominate the vacuum energy, and in order to avoid an 
unnatural hne tuning, the authors proposed [4I a subtraction of the Minkowskian vacuum energy 
based on the Hamiltonian formulation of General Relativity. 

Alternatively, it has been suggested that the use of a comoving cutoff could provide a way 
to renormalize the theory without introducing non covariant counterterms j^. Once more, the 
explicit calculation of the counterterms have not been worked out, and in fact, as we will see, the 
particular subtraction used in that work is not compatible with a redehniton of the bare constants 
of the theory. The comoving cutoff has also been considered to compute the effective potential for 
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B- 


interacting fields in curved backgrounds [6[1. In these works it is assumed that the Fourier modes 
that effectively contribute to the vacuum energy have a maximum frequency much lower than the 
UV cutoff that sets the validity of the theory. 


The use of an UV physical cutoff has also been advocated by other authors. In Ref.B h was 
assumed that the non conservation of the EMT is compensated by an additional non conserved 
source in Einstein equations. Alternatively, in Ref. B) the < Tqo > component was computed 
using a cutoff, while the other non trivial components < Ta > were determined from the SEE, 
in order to force the conservation of < > and therefore the consistency of the SEE. In these 

works, there is no discussion about subtraction: on the contrary, while the quartic divergences are 
cancelled due to the supersymmetry of the theory 


gravitate in Minkowski spacetime 


BJ, the quadratic divergences are taken as the physical values 


of the EMT. UV cutoffs have also been considered when performin g n umerical calculations in the 


context of non-equilibrium field theory in cosmological spacetimes [9|, and when computing loop 
corrections in inflationary models [r] |- 

In this work we will not consider other interesting alternatives, like the use of a covariant 


cutoff based on the Schwinger proper time approach 


lll |. or models in which the UV cutoff not 


only regulates the theory but also modifies the Lagrangian, giving rise to generalized dispersion 


relations for the quantum fields 


1 ^. 


The aim of this paper is to discuss in detail the renormalization procedure for quantum fields in 
cosmological spacetimes using as regulators comoving and physical three dimensional UV cutoffs. 
Perhaps not sufficiently stressed in the recent previous works on the subject BsL ™ addition to 
regularization, subtraction plays a crucial role in the renormalization process. Indeed, in order to 
absorb the infinities into the bare constants of the theory, the quantities to be subtracted should 
be generic, i.e., should not depend on the particular metric considered. We will show that, using 


adiabatic subtraction 


Q,B 


13l | it is indeed possible to introduce counterterms and end up with a 


consistent theory. However, the use of ultraviolet cutoffs make the choice of the counterterms rather 
unnatural. We will also show that, while quartic and quadratic divergences are strongly dependent 
on the regularization method, the logarithmic part is universal, and therefore the renormalization 
group flow of the gravitational constant is left unchanged. We will make a comparison of the 
resulting renormalized SEE, and show that they depend on the regularization scheme, and differ 
from the renormalized SEE obtained using the standard covariant procedure. 

The paper is organized as follows. In Section II we describe the SEE in the context of FLRW 
spacetimes. In Section III we describe the renormalization of the SEE using comoving and physical 
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cutoffs. We show that, for a comoving cutoff, it is necessary to introduce bare conserved fluids in 
order to absorb the divergences of the quantum EMT. For a physical cutoff, the divergences are a 
combination of non conserved tensors, and therefore it is necessary to include non covariant terms 
into the SEE to renormalize them. We discuss in detail the relation with the usual renormalization 
approach, and make some comments related with the cosmological constant problem. In Section 
IV we consider a Xcj)'^ theory, and discuss the renormalization of the equation for the mean value of 
the held, once more using non covariant methods. Section V contains the conclusions of our work. 

II. SEE IN ELRW SPACETIMES 

In a curved spacetime with a classical background metric the SEE read 

SttGb = {T^y) + (1) 

(1 2 ) 

where Gfj^y is the Einstein tensor, H^y and H^y are tensors that come from the variation of 
quadratic terms in the curvature in the gravitational action {B?, R^yR^^'^ and R^yp^jR^^^^^ , respec¬ 
tively), {Tpy) is the expectation value of the EMT of the quantum helds, and is the EMT of 
a perfect huid, that will be needed to renormalize the theory. The subindex B indicates that the 
corresponding gravitational constants, even those contained in are the bare ones. 

The expectation value (T^y) has UV divergences, that should be absorbed into the bare constants 
of the theory. In order to make this point explicit, we assume that (Tpy) is regularized in some 
way, and dehne its renormalized value as 

{Tpy) = {Tpy) — {Tpy) sub + (V) sub — (V) ren + {T^u) sub ; (2) 

where we have subtracted an appropriate tensor {Tpy)sub that cancels the divergences of {Tpy). 
Inserting Eq.([2|) into Eq.([T]) we obtain 

SttGb + (^SsHpy = {Tpy)ren + {Tpy)sub + ’ (^) 

and, after absorbing the divergences of {Tpy)sui, into the bare gravitational constants 

SttGr + Ol 3 RHpy = {Tpy)ren + ^{Tpy) + . (4) 

It is important to stress that one should include in the theory all the bare terms that are necessary to 
absorb the inhnities of {Tpy)sub- Moreover, depending on the choice of {Tpy)sub and of the regulator, 
it may happen that only a part of {Tpy)sub is absorbed into the bare terms. The remaining part 
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would give an additional contribution to the right hand side of the SEE, that we denoted by A{T^^) 
in Eq.dH). 

In the rest of the paper we will consider spatially flat FLRW background metrics, 

ds'^ = a{T)‘^{dT^ — dx"^) (5) 


where r is the conformal time. For these metrics, and in four spacetime dimensions, the tensors 
and H^i, are linearly dependent, and can be written in terms of [l|. Therefore, we will 
set Q !2 = 03 = 0 and denote oi = a. We stress that this is true only in four dimensions. Explicit 
expressions of the geometric tensors in FLRW metrics are presented in the Appendix. 

The EMT of a classical fluid has energy density and pressure given by 

p = Tlja^ p = T/Ja^ , (6) 


and the conservation equation reads 


p' + 3—{p + p) = 0 
a 

where a prime denotes derivative with respect to r. 

We will consider a free scalar field with classical action given by 


Smatter — 

SO the field equation reads 


I d‘x^( 




(p + rn^ + iR)(t) = 0 , 


( 7 ) 

( 8 ) 

(9) 


where ^ is the coupling to the Ricci scalar curvature. Defining (/> = y/a, the Fourier modes of the 
field satisfy 


Xfc + 


+ K “ 


a^R 


Xfc = 0, 


where uj\ = + m^a?. 

The EMT is given by 


Q 


= (1 - 2i)d^cl>du(t> + ( 2 ^ - 




( 10 ) 


( 11 ) 


The vacuum expectation value of the energy-momentum tensor can be written in terms of the 
Fourier modes of the field Xk- Taking into account the symmetries of the metric and vacuum 
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state, the non-vanishing components are (Too) and (T 
{Tii) = i((Too) — a^(T')j)). The explicit expressions are 


li = (T 22 ) = (T33) = {Tii). Therefore 




^^00) = J + WfcIXfcP + (^D{xkxi + xlXk) - 


{TV = 


= J dkfoo{k,T), 
1 


27r^a^ 


/ 


dk ¥ 


a^m^\Xk\^ ^ “ ^{XkX*k + Xlx'k) - ^^llxV 


a" o'2 


\Xk\^-^(C-l)(D' + lD^)\Xk\^^ = jdkf^{k,T), 


( 12 ) 


a a^j \ 6/\ 2 

where T = 2a' j a. The equations above have UV divergences coming from the high k behavior 
of the Fourier modes. In order to isolate the divergences, it is useful to consider the adiabatic 
approximation to the modes. Writing 

1 


Xkir) = 


Eg. dlOp is equivalent to the nonlinear equation 


exp — i 


Wk{v)dv , 


(13) 


T ..2 2 1 / 3 W'i 

^ 2 (Wk 2 Wi 


(14) 


In the large k limit, this equation can be solved iteratively using the number of derivatives of the 
scale factor as a small parameter. Doing this, and inserting the solution into Eq. (ll2|) . one obtains 
the adiabatic approximation to the EMT, {T^ad, that can be written as the sum of its divergent 
and convergent parts [l^ 

{TVad = (T^.)f)r + (WaT , (15) 


where 


/rT-\ \dlV _ 

{-^00/ad — 


87r2a2 / 


dk k^ 


1 


“^^k — —D I ~ 77 I I-1" 


1 




rn?a? 




U2 


(i^ 

16a;3 


-{72D"D - 36D'^ - 27D 


(T 


<a \div 
a/ad 


1 




/ 


dk k^ 


a?m? / 1 

" 6 


^k 


3D' 


+ 




UJ 




3D' + -D^ 
4 


li2 


(e-l) 


-(18D"' - 27D'D‘^) 


(16) 


The explicit expressions for the convergent parts are written in the Appendix. Defining the adia¬ 
batic order as the number of derivatives of the scale factor, by power counting it is easy to see that 
the zeroth adiabatic order diverges quartically, the second adiabatic order quadratically, and the 
fourth adiabatic order logarithmically. Higher adiabatic orders are finite. From these equations we 
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see that the adiabatic expansion is a powerful tool to isolate the divergences, whatever the regu¬ 
larization used. Moreover, as {T^u)ad can be written for any geometry in terms of time derivatives 
of the scale factor, it is the natural candidate for {T^u)sub- 


III. CUTOFF REGULARIZATION OF THE SEE 


A usual covariant approach to renormalize the theory in FLRW metrics is to combine dimen¬ 
sional regularization with adiabatic subtraction. Indeed, performing all calculations in an arbitrary 
number of dimensions n, one can show that the divergences in {T^y)ad are all proportional to the 
geometric tensors , Gfj^y , H^y and H^y . Therefore they can be absorbed into the gravitational 
constants Ab,Gb, and OiB- A subtle point in the renormalization process is that the fourth adia¬ 
batic order contains divergences that, as long as n 7 ^ 4, are proportional to geometric tensors and 
diverge as l/(n —4), renormalizing the coupling constants aBi, 2 - However, in the limit n —)• 4 they 
become hnite terms that modify {T^y)ren-, giving rise to the well known trace anomaly |l|. 

In this section we will explore the consequences of regularizing the theory imposing UV cut¬ 
offs in the momentum integrals. Our guiding principle will be that the divergences in {T^y)sub 
should be absorbed into appropriate bare terms in the SEE, and this implies that the subtraction 
should involve the adiabatic EMT up to the fourth adiabatic order. The last point needs further 
clarihcation. Using an UV cutoff, the renormalized EMT will be of the form 


(7; 


I'kuv 


pbV Iren 


dkT, 


I'kuv 


flU 


{k,T)- [ dkTp{k,T)=l dk{T^y{k,T)-Tp{k,T)), (17) 

Jo Jo 


where we made the reasonable assumption that {T^y)sub can be written as an integral over the 
momentum. The subtraction should cancel the UV divergences, and, at the same time, should be 
written in terms of the scale factor and its derivatives. The last property is crucial in order to 
absorb the divergences into bare constants of the theory. 

Note that when the subtraction is performed mode by mode, as in the last equality in Eq. (ll7p . 
it is possible to take the limit kuv 00 . 

As in the high k regime the modes of the field are well described by the WKB solutions to 
Eq. dlOp . the adiabatic modes will cancel the divergences. Therefore, the divergences of T^y will 
be cancelled by or at least by its divergent part. Moreover, in the adiabatic approximation 
the EMT becomes a local function of the scale factor and its derivatives, and this will allow for a 
redefinition of bare constants to cancel the divergences. 


A simpler procedure would be to subtract the contribution of the modes with kmax < k < kjjv 
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as proposed in Ref.^j. In this case, the renormalized energy-momentum tensor would be 

pkmax ^ 

{Tf,u)ren= / dk ff,^{k,T) . (18) 

JQ 

While this is obviously a finite quantity, the divergent part should be computed for each metric in 
terms of the modes that solve Eg. (1101) . and the result would not be a local function of the metric 
and its derivatives. Therefore, the divergences could not be absorbed into bare constants of the 
theory. 


A. Comoving Cutoff 


We will first consider a comoving cutoff {kuv = ^c)- After evaluating explicitly the integrals in 
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Ea. ffTHl) we obtain a4 

/rp \div _ 

\-^00/ad — 


+ 




m 


Ibvr^ 


-500X77^ log [Ac/m] 


3Ag 

Svr^o^ 


e- 


1 \ a 


/2 


327r2 


+ ^^16^2^ ^00 log[Ac/m] 


6/02 
( 1 ) 


+ 


87r2 


Goo log[Ac/m] 


m 


\div 
ii) ad 


1 At 


3 167r2o2 

A/ 


+ 


87r2o2 


1 A‘lm? 
3 167r2 


1 


m 


-5n^^log[Ac/m] 


6 




e-TT 2--3-^ + 


m 


Stt^ 


Gii log[Ac/m] 


(c _ l'\2 

+ , ^ o log[Ac/m], 


(19) 


167r2 

where we omitted finite terms. 

One can readily check that the adiabatic EMT is a linear combination of conserved tensors, each 
one with a different dependence on A^ and m. Eor instance, the zeroth adiabatic order corresponds 
to the sum of three fluids: a radiation fluid {p = p/3) proportional to A^, a fluid with equation of 
state p = —p/3 proportional to m^A^, and a cosmological constant {p = —p) with a logarithmic 
divergence. The quadratic divergence of the second adiabatic order corresponds to a fluid with 
energy density proportional to a'^/a®. The logarithmic divergences of all adiabatic orders are 
proportional to geometric tensors. 

The main points of this calculation are that the introduction of a comoving UV cutoff does 
not spoil the conservation of the EMT, and that it requires the introduction of non-standard bare 
fluids into the SEE in order to absorb the divergences. In what follows will compute explicitly the 


counterterms needed to renormalize the theory. The SEE with the appropriate bare fluid read 


1 


SttGb V 




+ AePoO + OisHo^ + /3B-n + 7s “ ^b3—^ = {Too)ren + {Too)ad 


1 




( 20 ) 















1 


8ttGb V ® 




+ ^BSii + + 


(i)^/3b1 1 


3 a2 


“ o7s + <Js 2^ - 3 


,'2 


9 

= {Tii)ren + {Tii)ad (21) 


As anticipated, we introduced three new bare constants in the theory. The relations between bare 
and renormalized constants read, using minimal subtraction 


A — A ^ 
Aren - Afi - ^ 


log[Ac/m] 


Gr 


1 


m 




vr 

n2 


log[Ac/m] 


(^ren — CHQ 
Pren — Pb 
^ren TB 
Sren — 


(i^ 

Ibvr^ 

A^ 


■log[Ac/m] 


Ibvr^ 

2A? 


m 


IGtt^ 

Stt^ 


( 22 ) 


The SEE without additional classical fluids can be obtained by setting j3ren = 0,jren = 0 and 

^ren — 0- 


B. Physical Cutoff 


Instead of using a comoving cutoff, one could alternatively use a physical cutoff kuv = o(r)Af/y, 
where Af/y is constant. The additional time dependence in the cutoff will compromise the con¬ 
servation of the adiabatic EMT, and therefore the renormalization of the SEE will involve non 
covariant counterterms. 

The divergent part of the adiabatic EMT is given by Eq. (jl9|) with the replacement Ac —)• aAjjv 

^uv'’ 


{Too) ad — 

+ 

+ 

{Tii)ad = 

+ 

+ 


Afr,c 


'■uv 


+ 


IGtt^ IGtt^ 


327r2 


log[Af/y/m] 


500 


A2 

Svr^ 




m‘ 


Stt^ 


log[A; 7 y/m] 


Goo 


log[Ac/y/m] 


1 Arn^rri^ 


1 ^UV I 

3167r2 3 


m 


167r2 327r2 


log[Ai7y/m] 


Qii 


A2 

87r2 




m 


87r2 


^^lG7r2^ log[Ac/y/m] , 


log[A; 7 y/m] 


A2 

47r2 




(23) 


where again we omitted hnite terms. 

The above equations clearly show the non-conservation of the adiabatic EMT. The quartic di¬ 
vergence corresponds to a fluid with p = l/3/o = constant, which is not conserved. The same 
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happens with the quadratic divergences of adiabatic orders zero and two. For example, the diver¬ 
gences of adiabatic order two consist of a conserved term proportional to the Einstein tensor 
plus a non conserved contribution proportional to Kuy5. Note however that all logarithmic 
divergences are proportional to conserved geometric tensors. 

In order to compute the counterterms, the SEE should be generalized introducing non-conserved 
bare terms (that can be thought as coming from the variation of non-covariant terms in the clas¬ 
sical gravitational action ^). Thus, in the ii-semiclassical equation we split the Newton and 
cosmological constants in two different constants, as follows 


+ Afig'oo + = (Too)ren + {Too)ad 


(24) 


1 


.o" a'2 


SttGb 

Inserting Eg. (1231) into Eqs 


a J SttGb « 


/2 


+ (A_B + ^B)gii + CUbH^- ^ — {Tii)ren + {Tu) 


ad 


(25) 


and (l2^ we obtain 


Afjy m^Kuv , 1 r. , . 


IGvr^ IGvr^ 


327r2 


Gr 


Gb 


TT 


\og[Kuv/m] 


^ren — 


(e-l) 


167r2 


A _ A 4 Af, 
^ren — “t“ 


■log[A; 7 V'/m'] 
2 


+ X- 


1 


Gr 


3167r2 3 167r2 

Gb 47r2 


(26) 


The covariance of the SEE is restored by setting to zero the renormalized values of the extra 
constants l/G^en = 0 and Aren = 0. 


C. Comparison with covariant regularization 


In the usual covariant approach (dimensional regularization plus adiabatic subtraction) one 
subtracts {Tfj_r,)sub = {T^y)ad-, where the adiabatic EMT contains all divergent and finite terms up 
to the fourth adiabatic order. The calculations should be entirely performed in n dimensions from 
the beginning, and in this case it can be shown that both the finite and divergent terms of {T^u)ad 
can be absorbed into the gravitational constants of the theory [l^. Moreover, when calculating 
{T^v)ren as the difference of two divergent integrals, the regulator can be removed performing the 
subtraction mode by mode. One of the consequences of this procedure is the appearance of the 
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trace anomaly, which is produced by the finite terms in {T^u)ad subtracted from (T^u)- Indeed, the 
trace of the renormalized EMT reads, for = 0 and ^ = 1/6 [],Q 




1 


/ n2'. 


OGOvr^a^ 


(D"' - D'D 


(27) 


while vanishes at the classical level. 

The situation in the presence of UV cutoff has some subtle points that deserve clarification. 
Up to now, we have not specified whether we subtracted the full adiabatic EMT or only its di¬ 
vergent part. In order to make contact with the covariant renormalization, we will subtract the 
full adiabatic tensor [l^. In this case, as the cutoff can be removed in the expression of {Tf^u)ren 
(see Eg. lfTTl) !. all regularization methods, covariant or not, give the same answer for {T^^)ren- This 
does not mean, however, that the SEE are all equivalent. Indeed, in the covariant approach, the 
complete {Tfj,i,)ad is absorbed into the gravitational bare constants. Therefore, {T^y)ren is the only 
quantum contribution to the SEE. When using UV cutoffs, only the divergent terms are cancelled 
by counterterms, and there are additional contributions coming from the finite parts of (T^j/)ad, 
that we denoted by A(T^i,) in Eq.Q. 

Eor a comoving cutoff we have 


A(roo) = 


4 2 

m a 
1287r2 


(1 -h41og(o/2)) - 


(e-i) 


m 


IGtt^ 


-(3 + 21og(a/2))Goo 


m 


1 


(c _ l')2 

-1^(1 + log(a/2))H^l^ - —^Goo- 

Iq^ 2 y ^ / jj 00 2887r2 153607r2o2 

it2 


+ 


1 


172807r2 00 

4„2 




2887r2 


00 ' 647r2a2 

2 / 


^Tu) = - 




3847r2 
it2 


(7 + 121og(a/2))- 


I67r2 
2 


(I -|- 2 log(a/2))Gj; 


^^^^^(l + log(a/2))i7|^) - "" 
vn?‘{^ — i) 9 


2887r2 


m2(^ — 1) 
ii ^ 


1 


327r2 

(e-l) 

2887r2 


57607r2a2 
2 


1 


D'D‘^ - + 


87r2 

1 


h: 


(1) 


^(1) _ (? - s) 


I927r2a2 


172807r2 “ 

45 


72D"D + 54H'2 -b 54T)'T)2- 

2 


(28) 


We stress that to compute A(T^j^) one should take into account not only the convergent part of 
the adiabatic EMT (see Appendix), but also the finite terms coming from One can check 

that A(r^iy) is a conserved tensor, and that the contributions proportional to log(a) are crucial for 
the conservation. 

Unless the finite extra terms described in Eq. (j28j) are artificially cancelled out by new coun¬ 
terterms, the resulting SEE differ from the usual ones, since the quantum effects are described by 
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the effective EMT = {T^u)ren + ^{T^y). For example, in the conformal limit (m^ = 0 and 

^ = 1/6), this effective EMT is traceless. Indeed, {T^y)ren has the anomalous trace Ea. dTFD . which 
is exactly cancelled by A(T/‘). This can be readily checked by an explicit calculation, but it is 
to be expected because, while the full {T^u)ad is absorbed into the bare constants when using a 
covariant regularization, this is not the case for an UV cutoff. 

The analysis can be repeated using a physical cutoff. The expression for A(T^j,) would be 
still given by Eq. (|28j) . omitting all terms proportional to log(a). Therefore, the effective energy- 
momentum tensor would not conserved. For consistency of the SEE, the non conserved terms must 
be cancelled by additional hnite counterterms. The remaining conserved terms read 

1 -r(l) ii-il rrll) 1 


A(roo) = - 

A(r,,) = - 


m 


2887r2 


m 


2887r2 


Goo + 


Gii + 


T-tinr\ 




172807r2 00 2887r2 oo I53607r2a2 

^ o) jj'(O) _j_ ^ 


172807r2 ** 


2887r2 ** 57607r2a2 




- -D^ 


(29) 


Note that the additional counterterms should cancel all terms proportional to — 1/6) and 

(^ — 1/6)^, since they produce non conserved contributions to A(T^,y) [l^. Note also that while the 
terms proportional to Gfj_y and Hj^J can be absorbed into a finite redefinition of the gravitational 
constants, the remaining contributions to A.{T^y) are nontrivial. The effective EMT for a physical 
cutoff does not have a trace anomaly and differs from the one computed using a comoving cutoff 
by local terms. 

D. A comment on the cosmological constant problem 

The zero point energy of the quantum fields gives an enormous contribution to the cosmological 
constant j^. A naive estimation in Minkowski spacetime consists in considering the sum of the 
ground state energy of each mode of the held. The sum is performed with a three dimensional 
cutoff to be of the order of Planck mass. Taking this naive estimation as the value of the cosmo¬ 


logical constant, it gives a disagreement o 


122 orders of magnitude with respect to the observed 


value. However, as pointed out in Ref. [l7|, the three dimensional cutoff does not respect Lorentz 
invariance, and the problem is alleviated when considering dimensional regularization, since this 
regularization kills the power law divergences, and keep only the logarithmic ones. Moreover, with 
a three dimensional cutoff in Minkowski spacetime, the vacuum expectation value of the EMT does 
not correspond to a cosmological constant but to that of a radiation huid. 

The results presented in this section show similar characteristics, generalized to curved space- 
times. When using a comoving cutoff, the renormalization of the cosmological constant involves 
only a logarithmic divergence. However, it is necessary to introduce into the theory three new bare 
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fluids, whose renormalized values should be fine tuned to zero. The quartic divergences usually 
advocated as contributions to the cosmological constant, now give naive estimations to the ampli¬ 
tude of the bare fluids. On the other hand, when using a physical cutoff, in addition to the usual 
quartic contribution to both cosmological constants Ar and A/j, it is necessary to fine tune the 
renormalized value of the new constants (A/j and 1/Gr) to zero in order to respect the covariance 
of the theory at the renormalized level. 

n 

These results reinforce the idea [18|] that the power law divergences should not be taken seriously 
as estimations of the zero point energy contribution to the cosmological constant. After absorbing 
the infinities into the bare constants, a finite piece coming from the logarithmic divergence remains, 
and is of order m^log[m/^], where ^ is an arbitrary scale (this is still enormous compared with 
the observed value of the cosmological constant, even for the electron mass). 


IV. MEAN VALUE EQUATION FOR SELFINTERACTING FIELDS 


We will now consider a scalar field with selfinteraction The classical action is 

Smatter = J > (^O) 


and the field equation 


0(j) [m\ + ^bR)4> + ~ ^ • 


(31) 


We have written the action in terms of the bare constants of the theory. 

We will be concerned with the renormalization of the equation for the mean value of the field 
00 = (0) • Writing 0 = 0o + 0, the equations for the mean value and the fluctuations of the field 
read tl9| 

(32) 


□00 + (m| + ^bR)4>o + ^00 + ^ 0 o( 0 ^) = 0 
□ 0 -|- ( mji + + “^000 = 0 - 


(33) 


The quantity (0^) is divergent and the infinities must be absorbed into the bare constants. Note 
that, as ( 0 ^) is already 0{h), we replaced the bare constants by the renormalized ones in the 
equation for the fluctuations and in the last term of the mean value equation. 

As for the EMT, in order to absorb the infinities into bare constants we will subtract the 
adiabatic expansion of (0^). We will insert the definition 

(0^)ren = (0^) “ (0^)ad 


(34) 
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into the field equation 

□ (/)0 + (w-B + ^bR)4>0 + ~^4>0 H- 7r4’o{(f^)ren H- TT 4’0 ad = 0 . 


,3 i /72\ I 1 / 12 \ 


(35) 


3 , T-u ' 2 """- ' 2 

and will analyze the divergences of {4>^)ad- For the same reasons as before, we choose the adiabatic 
expansion to perform the subtraction and renormalize {cjP'). A complete analysis using a covariant 


approach and for general metrics can be found in Ref. 


19l ]. Here we will restrict ourselves to the 


case of FLRW metrics, we will consider a time-dependent mean value (/>o(t) and will regularize the 
theory using UV cutoffs. 

The fluctuation held (/> is a free held with a variable mass = m\ + A/j(/>q/2 . As we are 
assuming that the mean value depends only on time, we can describe the held in terms its Fourier 
modes, that satisfy Eo. lfTOt) with uj‘l = + M'^a?. From Ea. (fT3]) we have 


47r2a2 / 


dk' 


Wk 


(36) 


The function Wk satishes Ea. (ll4p and can be solved using the adiabatic approximation. Although 
now not only the scale factor but also 4o depends on time, in principle one should include in the 
adiabatic expansion terms with derivatives of (po- However, it can be shown that these terms give 
hnite contributions to and therefore they can be omitted when discussing the renormaliza¬ 
tion of the mean value equation in the one loop approximation jl^ (there is no wave function 
renormalization for theory in the one loop approximation). 

Inserting the adiabatic approximation for Wk into Eg. (136]) we obtain [l^ 


(4- 


2\div 

ad 


1 


/ 


dkk^ 






(37) 


47r2a2 

In order to see explicitly the differences between the different regularization methods, we quote 
the result obtained within dimensional regularization. Replacing k'^ —>• ^”■“2 the integrand of 
Eq. (I37p and performing the integrations we obtain 

1 


22 \ div 
^ )ad 


87r2(n — 4) 




+ ~ 4 r 4 o + { Cr- 


R 


(38) 


where we omitted finite contributions. The hrst term renormalizes the bare mass niB, the second 
term the coupling constant A_b, and the third one the coupling to the curvature ^b- 
On the other hand, when using an UV cntoff kuv 


(4 


2\div 
/ ad 


1 

87r2 




1^1 log 


kuv 

afi 


(39) 


where /r is an arbitrary scale and we omitted hnite terms. 
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Let us first consider a physical cutoff kjjv = clAuv- From Eq. pQl) we see that the quadratic 
divergence renormalizes the mass, while the logarithmic divergence gives an additional term to 
the renormalization of the mass and renormalizes the other bare constants of the theory. The 
counterterms are, explicitly: 


= ^uv mllog{Auv/A) 

(40) 

^R = Cb ~ iog(A( 7 y//r) 

(41) 

3 

Afi = Ab - — A^log(Af/y///). 

Id 

(42) 


There are some similarities with the situation in the renormalization of the SEE. With an UV 
cutoff, the bare mass contains a quadratic divergence that is independent of the mass of the field. 
The logarithmic divergences has the same structure than in dimensional regularization. 

On the other hand, for a comoving cutoff we set kjjv = in Eq.([39|). We see that this choice 
presents an additional complication. The quadratic divergence depends on the scale factor, and 
therefore cannot be absorbed into a redefinition of the mass It is necessary to introduce an 
additional, non covariant term in the classical action of the interacting field 


^matter — 


(43) 


Assuming that the renormalized value of the new bare constant cjb is zero, that is, choosing 

ArA2 


aR = 0 = aB + 




(44) 


we obtain 


mR = ml- log(Ac//i) 

U = Cb- log(^c/M) 

3 

Ai? = Ab - —Xllog{Ac/fi) 


(45) 

(46) 

(47) 


which are similar to the counterterms obtained within dimensional regularization. 

It is interesting to remark that, as in the case of the SEE, the logarithmic divergences have 
always the same structure. As a consequence, the /3-functions of the theory will not depend on the 
regularization. On the contrary, the renormalized equation for the mean value (and therefore the 
effective potential for the scalar field) does depend on the regularization. Indeed, the analysis can 
be perfomed along the lines of Section IIIC, to check explicitly that the mean value equations in 
different regularizations will differ by local terms. 
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It is also worth to note that the presence of the term proportional io ub jo? in Ea. (H3]) introduce 
an ambiguity in the SEE, and could complicate its renormalization. Indeed, in order to obtain the 
SEE in FLRW spacetimes, one should write the action in terms of the scale factor and the lapse 
function N, and the extra factor l/a? can be written in different ways in terms of a and N. 


V. CONCLUSIONS 


In this paper we have analyzed the renormalization process for a scalar field FLRW sacetimes, 
using UV cutoffs. Our main findings are the following: 


• If one assumes that the infinities of the theory should be absorbed into bare constants of the 
theory, the subtraction is crucial and involves the adiabatic EMT. 

• The divergences that result in the case of a comoving cutoff are conserved tensors, and 
can be absorbed into the gravitational constants and additional bare fluids that should be 
introduced ad-hoc into the theory. 

• For a physical cutoff, the adiabatic EMT is not conserved, and it is necessary to introduce 
non covariant counterterms. Fine tuning of the finite parts of the counterterms is needed to 
restore covariance at the renormalized level. 

• All the regularizations give different answers for the effective EMT that appears in the SEE, 
even using the same subtraction. In particular, both the comoving and physical cutoffs do 
not produce an anomalous trace in the effective EMT. 

• The logarithmic divergences are always proportional to geometric tensors. 

• In contrast with what happens in the SEE for free fields, when considering the mean value 
equation for interacting fields, the use of a comoving cutoff needs the introduction of non 
covariant counterterms to renormalize the theory. 


Regarding previous works on the subject, we have showed explicitly that the renormalization 
suggested in Refs. [ 4 , 5] can indeed be pursued. However, we have described important aspects 
that have been overlooked before: on the one hand, adiabatic subtraction is crucial to absorb 
the infinities into the bare constants of the theory. One cannot make a subtraction based on each 


particular metric considered 6]. Moreover, even after renormalization, the final SEE are different 


than those coming from the usual approach, and different for comoving and physical cutoffs. On 
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the other hand, there is no reason to take the bare value of the EMT as its physical value, even 
after discarding the quartic divergences, as done in Refs. 7|, l8|. 

In summary, the use of UV cutoffs, although at first sight more physical than dimensional 
regularization, presents additional complications and arbitrariness, and in our opinion do not shed 
additional light on the cosmological constant and dark energy problems. 

We would like to end the paper with a comment related to the use of infrared (IR) cutoffs. 
Present knowledge establishes that the seeds for the inhomogeneities in the universe are quantum 
fluctuations that become classical when their wavelengths become larger than the horizon 211]. It 


is therefore natural to consider an effective held theory for the long wavelengths modes of the held, 
integrating out the short wavelengths. Therefore, the quantum part of the held is restricted to 
modes with k > Ha{t) and one has to face the problem of renormalizing the UV divergences in the 
presence of a time dependent IR cutoff. This is of course not restricted to cosmological applications. 
The same issue is present in the study of phase transitions, where the long wavelength part of the 
held associated to the order parameter of the transition becomes classical and non homogeneous, 
giving rise to the formation of domains ^ . A naive generalization of the usual adiabatic renormal¬ 
ization to the case of a time dependent IR cutoff gives of course a non-conserved energy-momentum 
tensor, even for free helds and using a covariant regularization. The non conservation is in this 
case physical, since during time evolution new modes enter into the classical held, and therefore 
there is energy exchange between the quantum and classical parts of the held. Only the full EMT, 
including the quantum, classical and stochastic parts should be conserved. Alternatively, one could 
consider comoving IR cutoffs 


2^. In both cases, dimensional regularization can be used to deal 


with the UV divergences. 


Appendix A 


In this appendix we collect some useful formulas used in the calculations. 
The explicit expressions for the geometric tensors in FLRW metrics are 

R=^(d' + ]-D 


Goo = 


Gii = D' 


~T 


r(l) 


9 /I 


= — 
a2 


,/2 


^00 = — 1 ^ ^ 


2D'" - D"D + -T)'2 - 3D'D^ + -D^ 
2 8 


(Al) 
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where D = 2a! ja. 

The convergent part of the adiabatic stress tensor mentioned in Ea. (|15p is given by [l 


Irp \conv _ 

\^00)ad — 

+ 


1 


Svr^a^ , 
7o®m® 
64a;9 


dkk^ 


a^rn^D'^ 


IGw® 64a;’^ 


[D'D^+D^)- -::^-- +(^-h 


1024a;ii 


{2D"D - T)'2 + AD'D"^ + T»^) 

(6T»"T» - 3T>'2 + 6T»'T)2) 


“8^ 


{T 


lO; \conv 
ex.) ad 


+ 


64a; 

1 

47r2a'^ 

32a;7 


{UOD'D'^ + 105T)'‘) + 




/ 


dkk^ 


8a;5 


(D' + D^) - 


64a;® 


+ 


(e-l) 


1 \2^2^2 


am 


(54T>'T)2 27T»'' 


8a;5 


32a;7 


(D'" + AD"D + 3D'^ + QD'D^ + D‘^) 

{2SD"D + 21T»'2 + UGD'D'^ + 49T>^) - m ^4^ 

25ocj^ 


128a;® 

1155ai0mi®T>4 / l\f9a^m^D‘^ ^ 9^,. 

+ —+ ( C - ^ ) ^^5 + ^(3^ + 6^ ^ + ^D'D^) 


2048a;i3 


6 


4a;'5 


4 2 

am 




„ (120T»"T> + 9011'^ + 390D'D^ + 105T>^) + „ 

32a;' 128a;® 


2 

i' 7 ~i2 


{IGSOD'D^ + 1365T>" 


945a®m®T>^ 


1\V 


n«-H 


4 4 

am 


128a;ii 

(270T>'T>2 + 135T)^ 


a?rin? 

32a;5 


(432T)"T) + 324D'2 + 648T>'T>'' + 27T> 


/ n2 


16a; 


(A2) 


ACKNOWLEDGEMENTS 


FDM would like to thank D. Lopez Nacir and L. Trombetta for interesting discussions on 
related matters, and ICTP for hospitality, where part of this work has been done. This research 
was supported by ANPCyT, CONICET, and UNCuyo. 


[1] N.D. Birrel and P.C.W. Davies, Quantum Fields in Curved Space (Cambridge University Press, Cam¬ 
bridge, England, 1982); S. M. Fulling, Aspects of Quantum Field Theory in Curved Spacetime (Cam¬ 
bridge University Press, Cambridge, England, 1989); L. Parker and D. Toms, Quantum Field Theory 
in Curved Spacetime (Cambridge University Press, Cambridge, England, 2009) 

[2] S. Weinberg, Rev. Mod. Phys. 61 , 1 (1989); S. M. Carroll, Living Rev. Rel. 4, 1 (2001). 

[3] L. Parker and S. A. Fulling, Phys. Rev. D 9, 341 (1974); S. A. Fulling and L. Parker, Annals of Physics, 
87 , 176 (1974). 

[4] M. Maggiore, Phys. Rev. D 83 , 063514, 2011; L. Hollenstein, M. Jaccard, M. Maggiore, and E. Mitsou, 
Phy. Rev. D 85 , 124031, 2012. 






























19 


[5] F. D. Albareti, J. A. R. Cembranos and A. L. Maroto, Phys. Rev. D 90, 123509 (2014). 

[6] A. L. Maroto and F. Prada, Phys. Rev. D 90, 123541 (2014). 

[7] N. Bilk, Phys. Rev. D 83, 105003 (2011). 

[8] D. Bernard and A. LeClair, Phys. Rev. D 87, 063010 (2013). 

[9] D. Boyanovsky, D. Cormier, H. J. de Vega, R. Holman and S. P. Kumar, Phys. Rev. D 57, 2166 (1998). 

[10] L. Senatore and M. Zaldarriaga, JHEP 1012, 008 (2010) 

[11] F. Sobreira, B. J. Ribeiro and 1. L. Shapiro, Phys. Lett. B 705, 273 (2011); see also N. Bilk, S. Domazet 
and B. Guberina, Phys. Lett. B 707, 221 (2012). 

[12] A. Kempf and J. C. Niemeyer, Phys. Rev. D 64, 103501 (2001); A. Kempf, Phys. Rev. D 63, 083514 

( 2001 ). 

[13] T. S. Bunch, J. Phys. A: Math. Gen 13, 1297 (1980). 

[14] This is not the case if one computes the function Wk in four dimensions and regularizes only the 
integration measure in momentum space. 

[15] Note that this does not mean loss of generality, since we are adding and subtracting in the 

SEE (see Eq.©. 

[16] As already stressed, for a physical cutoff there are no contributions proportional to log(a), that were 
crucial for the EMT conservation in Ea. (l28l) . 

[17] E.K. Akhmedov arXiv:hep-th/0204048 (2002). 

[18] J. Sola, Journal of Physics: Conference Series 453, 012015 (2013); J. Martin, Comptes Rendus 
Physique, 13, 566 (2012); C. P. Burgess, arXiv:1309.4133 (2013). 

[19] J. P. Paz and E. D. Mazzitelli, Phys. Rev. D 37, 2170 (1988). 

[20] This point has been noticed in Ref.jgj. 

[21] V. Mukhanov, Physical foundations of cosmology, Cambridge, UK: Univ. Pr. (2005); A. D. Linde, 
Particle physics and inflationary cosmology, Contemp. Concepts Phys. 5, 1 (1990). 

[22] P. C. Lombardo, F. D. Mazzitelli and R. J. Rivers, Nucl. Phys. B 672, 462 (2003). 

[23] F. Lombardo and F. D. Mazzitelli, Phys. Rev. D53, 2001 (1996); F. C. Lombardo and D. Lopez Nacir, 
Phys. Rev. D 72, 063506 (2005). 




